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ABSTRACT 

The fullerene Ceo has four infrared-active vibrational transitions at 7.0, 8.5, 
17.4 and 18.9 /im. We have previously observed emission features at 17.4 and 
18.9 /im in the reflection nebula NGC 7023 and demonstrated spatial correlations 
suggestive of a common origin. We now confirm our earlier identification of these 
features with Ceo by detecting a third emission feature at 7.04 ± 0.05 /im in NGC 
7023. We also report the detection of these three C 6 o features in the reflection 
nebula NGC 2023. Our spectroscopic mapping of NGC 7023 shows that the 18.9 
/im Ceo feature peaks on the central star and that the 16.4 /im emission feature 
due to polycyclic aromatic hydrocarbons peaks between the star and a nearby 
photodissociation front. The observed features in NGC 7023 are consistent with 
emission from UV-excited gas-phase C 6 o- We find that 0.1-0.6% of interstellar 
carbon is in Ceo; this abundance is consistent with those from previous upper 
limits and possible fullerene detections in the interstellar medium. This is the 
first firm detection of neutral Ceo in the interstellar medium. 

Subject headings: ISM: molecules — ISM: lines and bands — ISM: individual objects 
(NGC 7023, NGC 2023) — line: identification 
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Introduction 



Fullerenes are cage-like molecules (spheroidal or ellipsoidal) of pure carbon, such as 
C6Cb C 70 , C 76 , and C 84 . C 60 , also known as buckminsterfullerene, is the most stable fullerene 



and can account 



( jKroto et al. 



or up to 50% of the mass of fullerenes generated in the laboratory 



19851 ). Theorists have suggested that fullerenes might form around stars with 



carbon-rich atmospheres, su ch as carbon stars, Wolf-Rayet (WC) stars, and carbon-rich, 



hydrogen-poor R Cr B stars (|Kroto fc Jura 



2000 



Pascoli fc Polleux 



1992 



Goeres fc Sedlmayr 



1992 



Cherchneff et al. 



jon-rich grain 



2000). Fullerenes may also form as part of the car 
cond ensation process known to occur in the ejecta of Type II supernovae (jClayton et al 
20011 ). Hydrogenated amorphous carbon grains in the interstellar medium (ISM) may also 
decompose after interstellar shocks into polycyclic aromatic hydrocarbons (PAHs) and 



992 



fullerenes (IScott et al. 



Tielens 



( IBettens fc Herbst 



19971 ); PAHs comprise 9-18% of interstellar carbon ( jJoblin et al. 



2008). Ful 



1996 



erenes might also form via cold interstellar gas-phase chemistry 



19971 ). 



Foing &: Ehrenfreundl (119941 ) propose that two diffuse interstellar bands at 958 and 



963 nm are due to singly ion i zed C m , or C^ , altho ugh the identification is debated (IMaier 



1994 



Jenniskens et al. 



19971 ). 



Misawa et al. 



(120091 ) attribute additional diffuse interstellar 



bands at 902, 9 21, and 926 nm to C£ n - N 



post 



1995 



AGB stars (jSomerville fc Bellis 



Nuccitelli et al. 



Cami et al. 



1989 



o fullerenes were found towards carbon-rich 



Snow fc Seab 



20051 ). or R CrB stars (IClavton et al. 



1989), car 



1995 



Don stars ( 


Clavton et al. 


Lambert et al. 


2001 


)■ 



( 120101 ) have recently detected C 60 and C 70 in the planetary nebula Tc 1 (IC 



1266). 



Observational evidence for neutral fullerenes in the ISM, howeve r, has been e 



date. No ne utral fullerenes have yet b een found in the diffu se ISM (ISnow fc Seab 



Herbig 



usive to 



1989 



20001 ). dense molecular clouds ( iNuccitelli et al 



20051 ) . or at the photo dissociation 
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front in the reflection nebula NGC 7023 (IMoutou et al.lll999f ). 



Ceo has four infra red-active vibratio nal transitions, at 7.11, 8.55, 17.5, and 19.0 /im 



at 1200 K (gas-phase; 



Frum et al. 



(para-H2 matrix isolation; 



1991 



Sogoshi et al. 



. and at 6.98, 8.44, 17.3, and 18.9 /an at 2.4 K 



20001 ) . On this basis, we tentatively identified the 



. 7.4 and 18.9 /tm ISM emission features in the reflection nebula NGC 7023 as due to C 



( Werner et al. 



2004b 



Sellgren et al 



20071 ) 



60 



We provide here additional evidence for the presence of Cgo in reflection nebulae. We 
report here the detection of the predicted Ceo feature at 7.04 ± 0.05 /tm in NGC 7023. We 
also report the detection of Ceo features at 7.04, 17.4 and 18.9 /tm in a second reflection 
nebula, NGC 2023. The Ceo 8.5 /tm feature is too blended with strong 8. 6 /tm PAH 



emiss ion to be detected in reflection nebulae. We have found i n NGC 7023 (Sellgren et al. 



20071 ) that the 16.4 /tm emission feature, attributed to PAHs ( IMoutou et al 



2000|), has a 



spatial distribution distinct from that of the 18.9 /tm emission feature. We now compare 
the spatial distributions of the 16.4, 17.4, and 18.9 /tin emission features in NGC 7023, and 
find additional support for the Cgo identification. We also discuss the excitation mechanism 
for the infrared emission of Cgo- 



Observations 



We used the Spitzer Space 



Spectrograph (IRS; 



Houck et al. 



Tele scope ( j Werner et al. 



2004aJ) with the Infrared 



12004( 1 to obtain 5-38 /tm spectra of NGC 2023 (PI Sellgren, 
pid 40276, aorkeys = 23912704, 23911168) and NGC 7023 (PI Sellgren, pid 40276, aorkeys = 
23911424, 23911680). We obtained spectra with the short- wavelength low- resolution module 
SL (5-14 /tm; A/AA = 60-120), the long-wavelength low-resolution module LL (14-38 /tm; 
A/AA = 60-120), and the short-wavelength high-resolution module SH (9.5-19.5 /tm; A/AA 
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= 600). We chose nebular positions (29" west, 8" south of HD 37903 in NGC 2023; 25" 
east, 4" north of HD 200775 in NGC 7023) with a strong ratio of the 18.9 /im feature 
relative to the 16.4 /im PAH feature. We used matched aperture extraction in CUBISM 



(ISmith et al. 



2007bl ) to extract SL, LL, and SH spectra in regions of overlap between these 
spectral modules. The extraction aperture was 10.2"xl0.2" in NGC 2023 and 7.5"x9.2" in 
NGC 7023. 

We also retrieve from the Spitzer archive a spectral data cube for NGC 7023 with LL 
(PI Joblin, pid 3512; aorkey = 0011057920). We use CUBISM to derive spectral images in 
the 16.4, 17.4, and 18.9 /im features and 0-0 S(l) H 2 for NGC 7023. For the spectrum of 
each spatial pixel, we define a local continuum surrounding an emission feature or line, and 
subtract it before deriving the feature or line intensity. 

We search for bad pixels and correct them with CUBISM before extracting final 
spectra. We subtract dedicated sky spectra for the 5-38 /im spectra of NGC 2023 and NGC 
7023; no sky subtraction is done for the spectral mapping. 



Results 



Figured] illustrates our SL and SH spectra in NGC 2023 and NGC 7023. The 17.4 and 
18.9 /im emission features are prominent, and coincident with Ceo wavelengths. 

We show the 5-9 /im SL spectrum of NGC 7023 in Figure |2j We clearly detect an 
emission feature at 7.04 ± 0.05 /im. This feature is coincident, within the uncertainties, 
with the wavele ngth of the expecte d Cqq line. We highlight this emission feature by 



using PAHFIT (ISmith et al. 



2007aJ) to fit the 5-9 /im spectrum with a blend of PAH 
emission features in addition to the new emission feature at 7.04 /im. The full-width at 
half-maximum of the 7.04 /im Ceo feature is 0.096 ± 0.012 /im, significantly broader than 
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Wavelength (|a,m) 



Fig. 1— Spttzer-IRS spectra (solid curves) of NGC 7023 (25" east, 4" north of HD 
200775; top ) and NGC 2023 (29" west, 8" south of HD 37903; bottom), obtained with 
the short-wavelength low resolution module (SL; 5.2-10.0 /xm; A/AA = 60-120) and the 
short- wavelength high-resolution module (SH; 10.0-19.5 /iin; A/AA = 600). We mark Ceo 
lines at 7.04, 8.5, 17.4 and 18.9 /im (vertical lines). The strong emission feature at 8.6 /im 
is due to PAHs. H 2 emission lines fall at 9.66, 12.3, and 17.0 /xm. 
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our spectral resolution. We also detect the 7.04 /im Ceo feature in NGC 2023. We present 
the C 60 band intensities in Table [TJ 



In our previous long-slit spectroscopic investigation of NGC 7023 (ISellgren et al.ll2007l ). 
we found that the 18.9 /im feature peaks closer to the central star than PAHs. We now 
illustrate this more clearly with the LL spectroscopic map extracted in NGC 7023 (Figure 
[3]). The 18.9 /im emission is clearly centered on the star. By contrast, the 16.4 /mi PAH 
emission peaks outside the region of maximum 18.9 /im emission, in a layer between the 
star and the molecular cloud. The photodissociation front at the UV-illuminated front 
surface of the molecular cloud is delineated by 0-0 S(l) H 2 emission at 17.0 /im. 



Our previous observations (ISellgren et al.l 120071 ) suggested that the 17.4 /im feature 
might be a blend of a PAH feature and an emission feature with the same spatial distribution 
as the 18.9 /im feature. We now confirm that this is the case with IRS/LL spectroscopic 
imaging. We show an image of the 17.4 /im emission from NGC 7023 in Figure HI overlaid 
with contours of 18.9 /zm and 16.4 /xm emission. The 17.4 /im emission clearly shows one 
peak on the central star, coincident with 18.9 /xm Cgo emission, and a second peak co-spatial 
with 16.4 /im PAH emission. Thus there is an ISM component with emission features at 
17.4 and 18.9 /im, which has a different spatial distribution than PAHs traced by the 16.4 
/im feature. 

Our imaging spectroscopy demonstrates the spati al separation be t ween regions of peak 



PAH emission and peak Ceo emission (Figs. [3] and HI). 



Boersma et al. 



(120101 ) find that the 



16.4 /im feature f rom different sources corre 
feature does not. 



ates with other PAH features but the 18.9 /im 
Velusamy fc Langerj (j2008l ) observe the 18.9 /tm feature to peak on the 
central star in the reflection nebula NGC 2316 (Parsamian 18). They also find that the 17.4 
and 18.9 /im features have distinct spatial distributions in this object, and argue from this 
that Ceo is unlikely to be the carrier of these features. Our result that the 17.4 /im feature 
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Fig. 2. — Spitzer-IRS 5-9 /im spectrum of NGC 7023 (open squares), obtained with the 
short- wavelength low- resolution module (SL; A/ A A = 60-120). We mark Ceo lines at 7.04 
and 8.5 /im (vertical lines). We show the individual contributions of PAH features at 5.3, 
5.7, 6.2, 6.4, 6.7 , 7.4, 7.6, 7.8, 8.3, a nd 8.6 /im to the spectrum, by decomposing the spectrum 



with PAHFIT (ISmith et al. 



2007a]) and then overplotting the Drude profile of each feature 
(blue curves). The Drude fit to the Ceo feature we detect at 7.04 ± 0.05 /im is highlighted 
(magenta curve). The 8.5 /im Ceo feature is blended with the strong 8.6 /im PAH feature. 
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Table 1. Observed 3, and Calculated 13 Cgo Intensity Ratios 



-^7.04/ 1 



18.9 



hx/I: 



18.9 



Ivta/I 



18.9 



Object 

NGC 7023 (A/AA=60-130) 0.82 ± 0.12 
NGC 7023 (A/AA=600) 
NGC 2023 (A/AA=60-130) 0.086 ± 0.004 
Absorbed photon energy 

5 eV 
10 eV 
15 eV 



0.46-0.58 0.41-0.43 
0.76-0.94 0.57-0.59 
0.97-1.20 0.67-0.71 



0.42 ± 0.02 
0.33 ± 0.01 
0.47 ± 0.01 

0.28-0.38 
0.28-0.38 
0.29-0.38 



a Observed intensity ratios, derived using PAHFIT (jSmith et al.ll2007ai ). 
We give statistical uncertainties; systematic fitting uncertainties are 15% 
for the 7.04 /xm intensity ratio and 30% for the 17.4 /jm intensity ratio. 
The observed 17.4 fim feature has not been corrected for PAH emission 
blended with it. 



b E mission spectrum calculated with Monte Carlo code ( jJoblin et al. 
20021 ) for molecular co oling cascade after absorb i ng a stellar photon . C 60 



vibrational data from 



Schettino et al. 



(2001) 



Menendez k Page! fl2000h . 



Choi et al. 



(12000h . and 



^Intensities (W m 2 sr 1 ) normalized to the 18.9 /im feature intensity 




Fig. 3. — Three-color image of NGC 7023, in 18.9 /zm Cgo emission (red), 16.4 fim PAH 
emission (green), and 17.0 fim 0-0 S(l) H 2 emission (blu e). We constructed each image by 



2007bl ). We illustrate 



analyzing Spitzer /IRS-LL long-slit spectra with CUBISM (ISmith et al.l 
where we measured the 5-38 fim spectra shown in Figs. [1] and [2] ( white rectangle ). We mark 
the location of HD 200775 (star). Each pixel is 5.1" x 5.1". The 18.9 fim Cqq feature 
peaks on the central star, while the 16.4 /zm PAH emission is brightest between the 18.9 /zm 
emission region and the photodissociation front traced by H 2 emission. 



17.4 




Fig. 4.— Image of 17.4 /im emission in NGC 7023 (g ray-scale image), d erived from 



Spitzer /IRS-LL long-slit spectra analyzed with CUBISM ( pmith et al. 



2007bf ). We over- 



lay contours of 18.9 /im Ceo emission (red) and 16.4 /im PAH emission (green). Each pixel is 
5.1" x 5.1". One component of the 17.4 /im emission is co-spatial with 18.9 /xm Cgo emission 
and the other component of the 17.4 /im emission follows 16.4 /im PAH emission. This 
illustrates that the 17.4 /im emission feature is a blend of a PAH feature and 17.4 /im Cgo 
emission. 
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is a blend of a PAH feature and Ceo i n NGC 7023 naturally explains their observations. 



4. Discussion 



To derive the abundance of Ceo, we assume that fullerenes and PAHs absorb UV 
starlight, and then re-radiate it all as infrared emission features. We do not include any 
potential visible fluorescence by either Ceo or PAHs. For fulleren es, we adopted experime ntal 



values of the ab sorption cross-section of C 6 o at 0.09-0. 35 um (lYasumatsu et al. 



Yagi et al. 



Yagi et al. 



2009 



. The absorption cross-sections from 



Yasumatsu et al.l ( 



1996 



19961 ) and 



(120091 ) differ by a factor of 2, indicating the overall uncertainty in our abundance 



calcu 
from 



ation. For PAHs, w e adopted the absolute absorption cross-section at 0.09-0.30 /im 



Li fc Draind (120011 ). We integrated both of these absolute absorption cross-sections 



over a 17,000-22,000 K blackbody, as is appropriate for the central stars of NGC 7023 and 
NGC 2023. We find that Ceo an d PAHs have similar integrated UV absorption strengths 
per C atom. 

We compare the sum of the intensities of the 7.04, 17.4, and 18.9 /xm features, assumed 
to be due to Ceo, to the sum of all other infrared emission features at 5-38 um, assum ed 



2007al ) to 



to be due to PAHs. We analyze SL and LL spectra with PAHFIT (jSmith et al. 
find that the ratio of Ceo to PAH emission is 0.01-0.03 in our observed positions in NGC 
7023 and NGC 2023. By adopting a percentage of interstellar carbon in PAHs of 9-18%, 
we derive a percentage of interstellar carbon C 6 o, p(Ceo), of 0.1-0.6% in regions of bright 
Ceo emission. 

Our p(Cqo) value i s consistent with othe r estim ates of p(Ceo) and of the percentage of 



carbon in Cg , p(Cq ). iFoing fc Ehrenfreundl (119941 ) estimate p(Cq ) = 0.3 -0.9% from tw o 
diffuse interstellar bands at 958 and 963 nm which they attribute to C^~ . iHerbig! (120001 ) 
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use these same two bands to estimate p(Cqo) + = 0.1-0.3% in diffuse clouds. iHerbigl ( 120001 ) 
places an upper limit of p{C m ) < 0.0008% in diffuse clouds, showing that C 60 is primarily 



ionized in diffuse clouds. 



Moutouetal 



(119991 ) measure p{C m ) < 0.3% and p(C m ) + < 0.3% 



from the lack of emissio n features at 7.0-8.5 /i m, at a position in NGC 7023 where the 18.9 



/im Cqq feature is weak. 



Nuccitelli et al. 



(120051 ) find p(Cqo) < 0.6%, from their non-detection 



of the 8.5 /im feature in absorpti on towards R CrB (HR 5880) and three massive young 



stellar objects. 



Cami et al 



( I2010l ) derive p(Ceo) = 1.5% in planetary nebula Tc 1. 



The relative intensities of the Cgo bands provide information on the conditions in which 
the molecule emits. To probe these conditions, we use a Monte Carlo code, based on a 



micro-canonical formalism, devel oped to simulate 



the absorption of a UV photon (jJoblin et al. 



; he emission cascad e of PAHs following 



2002 



Mulas et al. 



20061 ). We calculate the 



evolution of the internal energy in the molecule and the number of photons emitted in 
each infrared band during the cooling cascade. We calculate the emission spectrum after 
C 60 absorbs a 5, 10 or 15 eV photon, reaching temperatures of 8 00, 1200 or 1570 K at the 
beginning of the cooling cascade. We use the list of modes from iMenendez fc Paed (120 00') 



and infrared intensities fat K) calculated using density functional theory (IChoi et al 



2000 



Schettino et al. 



200l|). 



We compare the predicted line intensities with the observations in Table [TJ We find 
the ratio of the intensities of the 17.4 and 18.9 /im bands is not sensitive to the energy of 
the absorbed UV photon. The observed ratio varies but this is likely because the 17.4 /im 
band is blended with a PAH feature. The ratio of the intensities of the 7.04 and 18.9 /im 
bands, however, is very sensitive to the absorbed UV photon energy. Table [1] shows the 
7.04 /im intensity in NGC 7023 is consistent with the cooling cascade of Ceo excited by UV 
photons with a mean energy of 10 eV. A lower photon energy is suggested for NGC 2023, 
perhaps because of its blister geometry with the star in front of a dense molecular cloud. 
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If Ceo in NGC 2023 is excited by lower energy photons, then this would also affect the Ceo 
abundance derived there. More detailed modeling will be needed to clarify this. 

Our interpret ation of the NGC 7023 Ceo emission differs from the very recent results of 



Cami et al. 



(120101 ) who concluded that Ceo in Tc 1 emits in solid phase at a temperature of 



330 K. Their analysis, however, is based on an excitation diagram to derive a single emission 
temperature, which is not appropriate to describe the cooling cascade of UV-excited 
molecules. 



5. Conclusions 

We c onfirm our identificat i on of the 17.4 and 1 8.9 fim emission features in NGC 7023 
2004bl 



with Ceo (I Werner et al. 



Sellgren et al. 



20071 ) by detecting a predicted emission 



feature at 7.04 ± 0.05 /im. We also detect 7.04, 17.4, and 18.9 /zm emission features in 



NGC 2023. We demonstrate that the 17.4 /im emission feature in N GC 7023 is a b 



Cf;n and PAH emiss ion with Spitzer imaging spectroscopy. Our work (I Werner et al. 



end o 



2004b 



Sellgren et al. 



20071 ; this paper) is the first firm detection of neutral Cqq in interstellar space. 



We find that the infrared emission in NGC 7023 is consistent with the emission 
of gas-phase Cgo excited by UV photons. Further modeling is required to explain the 
observations of NGC 2023. The percentage of interstellar carbon in Ceo is 0.1-0.6% in NGC 
7023. This is consistent with previous estimates of, and limits on, the interstellar C 60 and 



Cg abundances. 



We thank Nick Abel, Lou Allamandola, Bruce Draine, Alain Leger, and Farid Salama 
for useful conversations, Dominique Toublanc for support with the Monte Carlo code, and 
Mike Jura for suggesting the Cqq identification. This work is based on observations made 
with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, 
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